Hyperon bulk viscosity in strong magnetic fields 



Monika Sinha and Debades Bandyopadhyay 

Theory Division and Centre for Astroparticle Physics, 
Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata-700064, India 

Abstract 

We study the bulk viscosity of neutron star matter including A hyperons in the presence of 
quantizing magnetic fields. Relaxation time and bulk viscosity due to both the non-leptonic weak 
process involving A hyperons and direct Urea processes are calculated here. In the presence of 
a strong magnetic field of 10^'^ G, the hyperon bulk viscosity coefficient is reduced whereas bulk 
viscosity coefficients due to direct Urea processes are enhanced compared with their field free cases 
when many Landau levels are populated by protons, electrons and muons. 

PACS numbers: 97.60.Jd, 26.60.-c, 04.40.Dg 
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I. INTRODUCTION 



R-mode instability plays an important role in regulating spins of newly born neutron 
stars as well as old and accreting neutron stars in low mass x-ray binaries jl|. Gravitational 
radiation drives the r-mode unstable due to Chandrasekhar-Friedman-Scliutz mechanism 
, llO| • R-mode instability could be a promising source of gravitational 
radiation. It would be possible to probe neutron star interior if it is detected by gravity 
wave detectors. 

Like gravitational radiation, electromagnetic radiation also drives the r-mode unstable 
through Chandrasekhar-Friedman-Schutz mechanism. There exists a class of neutron stars 



called magnet ars 



with strong surface magnetic fields 10^"^ — 10^^ G as predicted by 



observations on soft gamma-ray repeaters and anomalous x-ray pulsars [l2|, ll3|] . The effects 
of magnetic fields on the spin evolution and r-modes in protomagnetars were investigated by 



different groups 



14 
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161]. On the one hand, it was shown that the growth of the r-mode 



due to electromagnetic and Alfven wave emission for strong magnetic field and slow rotation 
could compete with that of gravitational radiation |l^. On the other hand, it was argued 
that the distortion of magnetic fields in neutron stars due to r-modes might damp the mode 
when the field is ~ 10^^ G or more Q, Q. 



The evolution of r-modes proceeds through three steps [I7j]. In the first phase, the 
mode amplitude grows exponentially with time. In the next stage, the mode saturates 
due to nonlinear effects. In this case viscosity becomes important. Finally, viscous forces 
dominate over gravitational radiation driven instability and damp the r-mode. This shows 
that viscosity plays an important role on the evolution of r-mode. Bulk and shear viscosities 



were extensively investigated in connection with the 
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23, l28|, [29, 
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damping o 
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34 



;he r-mode instability 



35l . l36l |. In particular, it 



was shown that the hyperon bulk viscosity might effectively damp the r-mode instability 
251]. However all these calculations of viscosity were performed in the absence of magnetic 
fields. The only calculation of bulk viscosity due to Urea process in magnetised neutron 
star matter was presented in Ref . 37|] . This motivates us to investigate bulk viscosity due to 
non-leptonic process involving hyperons in the presence of strong magnetic fields. It is to 
be noted that the magnetic field in neutron star interior might be higher by several orders 
of magnitude than the surface magnetic field |38|]. Further it was shown that neutron stars 
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could sustain strong interior magnetic field ~ 10^^ G |39l. |40[|. 

The paper is organised in the following way. In Section [TTl we describe hyperon matter in 
strong magnetic fields. We calculate bulk viscosity due to the non-leptonic process involving 
A hyperons and due to leptonic processes in Section lllli We discuss results in Section IIVI 
and a summary is given in Section |Vl 



II. HYPERON MATTER IN MAGNETIC FIELD 



We describe 13 equilibrated and charge neutral neutron star matter made of neutrons, 
protons, A hyperons, electrons and muons within a relativistic mean field approach 
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The baryon-baryon interaction is mediated by a, and p mesons. In the absence of magnetic 



field, the baryon-baryon interaction is given by the Lagrangian density 
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44| 



B=n,p,k 



+ -{d^ad^a-mla^)-U{a) 



1 



-^P,. ■ P'" + 2^>M ■ 



The scalar self interaction term 

U{a) = 

and 



43 



44 



45| is. 



^ 9i rriN {goNCff + ^92 {gaNf^Y 



(1) 



(2) 



Pi,u = d^Pi, - d^p^ . 
In mean field approximation, the effective mass of baryons B is 

m*j^ = niB - gcrBO- , 
where a is given by its ground state expectation value 



0- = -\ ly^gcrB ris 



dU 



(3) 
(4) 

(5) 



(6) 
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The scalar density is given by 



n 



B 



m 



B 



B 



(2vr)3 7o 

The chemical potential for baryons B is 

l^B = \fkj~+m^ + 9u^B + pI QpB hB, 
where I^b is the isospin projection and 

n 1 \^ 

2^9u>B TLB, 



(7) 



—^y^SpB hB riB. 



(8) 

(9) 
(10) 



P B 



The total baryon number density is Ub = XIb'^b- 

Now we consider the effects of strong magnetic fields on hyperon matter. The motion 
of charged particles in a magnetic field is Landau quantized in the plane perpendicular to 
the direction of the field. We solve Dirac equations for charged particles using the gauge 
corresponding to the constant magnetic field Bm along the z axis as Aq = 0, A = (0, xB^, 0). 
In the presence of a constant magnetic field, the Lagrangian density for protons is taken 



from Ref. 



461 . The positive energy solutions for protons are 



1/2 



with ^ = \/b (^x — j and b = qBm- 
The positive energy spinors, lAy{k,x) 



^-i{et - kyy - k.z) U^^^^k^x)^ 



(11) 



47 



481 . 149l . 1501 1 are given by 



U^^,y{k,x) = 



e + m* 







e +m* "v^^ 



\ e'+m* 



2ub Tj 



u+l 



(0 



(12) 



and 



Ui^^{k,x) 



e +m* 







(13) 



where e' = y^p^ _|_ _|_ 2uqBm- 



The proton number density Up and scalar density are given by 46 1 

qB, 



27r2 

i/=0 



^^^fcp(z/), (14) 



= 4r-7r'TlI; In — ^ , (15) 



where = y kp^ + m*^ and kp{h') = ■sjkp^^ — 2uqBm- Maximum number of Landau 
levels populated is denoted by Vmax and the Landau level degeneracy g^, is 1 for z/ = and 
2 for > 0. Similarly, we treat noninteracting electrons and muons in constant magnetic 
fields. 

The total energy density of neutron star matter is 



B=nA ^ ^ 



9u kp{u)Hp + (rrip + 2uqBm) In ■ 



Similarly the total pressure of the system is given by 

2Jb + 1 f''^^ k^ dk 

B=n,A 



1 2 2 rr/ N 1 2 0^ 1 2 0^ 1 ^ 2Jb + I f''^^ k^ dk 



+ ] ^p('^)/^P ~ ("^P +2z/g5^)ln. 



where A;i(z/) = ykp^ — 2vqBrn . The relation between pressure and energy density defines 
the equation of state (EoS). 
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III. BULK VISCOSITY 



The macroscopic compression (or expansion) of a fluid element leads to departure from 
chemical equilibrium. Non-equilibrium processes cause dissipation of energy which is the 
origin of bulk viscosity in neutron stars. Weak interaction processes bring the system back 
to equilibrium. In this calculation, we consider the non-leptonic reaction 

n + p — > p + A , (18) 

as well as direct Urea (dUrca) processes which are represented by 

n — ^ p + r + ui, (19) 

where I stands for e or fi. When the chemical equilibrium is achieved, chemical potentials 
involved in above reactions satisfy yUn — /^a = and fin — f'p — fJ'i = respectively. In this case 
the forward and reverse reaction rates, Tf and F,. are same. The departure from chemical 
equilibrium due to macroscopic perturbation gives rise to the difference between forward and 
reverse reaction rates, F = — 7^ 0. For a rotating neutron star, the r-mode oscillation 
provides the macroscopic perturbation which drives the system out of chemical equilibrium. 



The real part of bulk viscosity coefficient can be written as 51 1 



^ l + (cur)2 \dnr,) duk ^ ' 

where Xi = Ui/nb is the equilibrium fraction of i-th species, u is the angular velocity of (/, m) 
r-mode and r is the microscopic relaxation time. For a neutron star rotating with angular 
velocity Q, the angular velocity {u) of (/, m) r-mode is given by 

2m 



UJ 



-Q. (21) 



/(/ 

We are interested in / = m = 2 r-mode in this calculation. The relaxation time is given by 

r 5fxnbSxn ^ 
where 6fi refers to the chemical imbalance. Here F is the total reaction rate. 

The partial derivative of pressure with respect to neutron number density can be evaluated 
from the EoS under consideration as 

op 7,2 gg-iv "^71 ^ 



The total derivative dxn/dri}, can be evaluated numerically. 

Now, we calculate relaxation times for above mentioned processes in presence of magnetic 
field Bm using the EoS as described in section [Tll 



A. Non-leptonic process 

Here we consider the non-leptonic process given by Eq. ( fT8l) . In this case, only protons 
are affected by magnetic fields. The reaction rate is given by 



Vd^kr, f L^dk„. Z""'' Lydk^,,^, f Lzdk„, /"''a'' Lydk^,, f Vd^k 



(27r)3 J 2tt J_m^ 2tt J 2n J_bi^ 2tt J (27r)3 

X F(e„,ep,,ep^,eA), (25) 

kp.^ and kp^^ being the z component of momenta of initial and final protons respectively and 
kn and /ea denote momenta of neutrons and A hyperons. The Pauli blocking factor is given 
by 

F(e.,e,,,ep„eA) = /(e.)/(e,J{l-/(e,,)}{l-/(eA)} - /(eA)/(e,,){l - /(e,J}{l - /(e.)} , 

(26) 

with the Fermi distribution function at temperature T 

fie.) = (27) 

1 + e kT 

The matrix element Wfi is given by 

Wf. = ^.^,1 , \M\' e-«^ 5{e)6{ky)S{k,) , (28) 

where 

= ~ ~ ^"''^^ and 6ik) = 6iK, + kp^-kp^-kA). (29) 
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The invariant amplitude squared for the process is 

+ {K ■ kp^)iKf ■ ^a){(1 + ^np)(l + ^pa) + %np^pA} 

+ {kn ■ kA){kp^ ■kp^.){{l + glp){l + g^^) - Ag^pgpA}] . (30) 

In calculating the matrix element given by Eq. fl28|) we use the solutions of Dirac equation 
for protons in magnetic field given by Eqs. (fT2|l and (fT3|) . We also assume that the magnetic 
field is so strong that only zeroth Landau level is populated. Now we integrate over kp^^ and 
kp^^ using S{ky) and obtain 

(27r)7V16 'J "7 '^^J ^\e^ep^ep^.eA),^,^^ 

y,^-[(k..^k,.)Hik„y~k,,)-]/2l> p^^^^ ^^^^ ^^^^ ^(^^ ^ _ _ ,^)5^k,)_ (31) 

Here the subscript 5{ky) denotes that this condition has been imposed on the quan- 
tity within the parenthesis. Next we perform the integration over fc„ and kp^ and write 
d^k = k"^ dk d{cos6) d(j). The delta function of z-components of momenta implies 
knz + kp.^ = kp^^ + kAz- Here we note that when protons occupy only the zeroth Landau 
level, they have momenta along the direction of magnetic field i.e. in z direction. Hence we 
have kpz = kp^. Then depending upon whether the initial and final protons are moving in 
the same or opposite direction we have /eaz ~ knz = or kp^^ ~ knz = '^kp^. Next we perform 
the angle integrations using 5{kz) and change variable /c to e to get 

b f k 



^F{en, tp^.tpj, eA)(5(e„ + tp^ - tp^ - ca). (32) 



Here the subscript 9int denotes the angle integrated value. As particles reside near their 
Fermi surfaces in a degenerate matter we replace momenta and energies under integration 
by their respective values at their Fermi surfaces. 
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The matrix element squared is rewritten as, 

ii\M\')o.Jsiky),5ik.) = ^Glsm'2e,[2m:m;'ml{l-gl^)il-glj,) 
- m;m;/ip/iA(l - glp){l + ^pa) - m;m>„/ip(l + glp){l - gl^) 

+ /in/^J/iA{(l + 9lp){l + qIa) + 4:9np9pA} 



Pi "^Pf 



{kTf 



As 6fi << kT, the energy integration of Eq. (!32l) can be written as 511] 

J dendep^dep^deA F{en, Cp^, e^^, eA) S{en + e 
Finally we get 

1 qBmkpj, 



27r2 



(33) 



(34) 



3847r3 fcl 



{i\M\')eJ 



2 /1l2 



e{(A; 



-pF„-fcFj2]/26 



+e{{kF^ - kp.y - Akj^je 



(A;T)2 (35) 



The expression of the reaction rate for a zero magnetic field is given by 

1 



5ll 



1927r= 



(1^1^) kp, {kry 



(36) 



where the angle averaged matrix element squared is same as given by 511] • 

Now the quantity 6fi/6xn in Eq. (12^ is to be evaluated under the condition of total 
baryon number conservation 51] 



Sun + Sua 



0, 



which leads to 



Further we have 



TT 



- CtnA - "An + "AA, with tti 



nia J fl* 1 i 



(37) 



(38) 



D 



H T9uji9u,j H J 5'pi^3i 5'pi^3i- (39) 



m 



m 



Here is the same as given by Eq. (1241) . Next we evaluate the relaxation time of the 
non-leptonic reaction at a given baryon density using Eq. fl22]) along with Eqs. fl5Sl) . fl55]) 
and (15^. 

As soon as we know the relaxation time, we can calculate the bulk viscosity coefficient C 
due to the non-leptonic interaction at a given baryon density from Eq. (!20|) . 



B. Leptonic processes 



Here we consider dUrca processes involving nucleons, electrons and muons in a magnetic 



field. The forward reaction rate for this process is then given by 



43,148 



49| 



V(fkn f VSky f L^dk^p f 2" Lydkyp f L^dk^i f 2" Lydkyi 



(27r)3 J (27r) 
xF(e„,ep,e/). 

Here F(e„, e^, ei) is given by 



27r 



Mi 27r 

2 



27T 



27r 

2 



F(e„,e„Q) = /(e„){l-/(ep)}{l-/(Q)}. 



(40) 



(41) 



Using the solutions of Dirac equations for protons and electrons in magnetic field, we obtain 
the matrix element 



W 



fi 



(27r)= 



■\M\^5{e)5{ky)5{k, 



V^LyL^) 

Firstly we treat the case following the prescription of Baiko and Yakovlev 



(42) 



48l | when 



protons and electrons populate large numbers of Landau levels. In this case, we have 

^\M\^ = 2Glcos'ec{l + ?>G\)F\ (43) 



where F is Laguerre functions for both protons and electrons 
rate is given by, 

where 



2>2tiGI cos2 e, 



(2vr) 



48| . The forward reaction 
'^"'"^^^'i?-|rf6.e^J(6.), (44) 



R"^ = 2 




1 

d cos 6pd cos 61 — ^ — -Fjj^^N^{u)Q{kF^ - {kp^ cos 6p + kp^ cos 6i\) , (45) 



-1 



and 



dendepdeiF{en, Cp, ei)(5(e„, - tp - ei - 
{kTf + {eJkTf 



(46) 



2 1 + e'-l^'^ 

As there is chemical imbalance due to the perturbation, the reverse reaction rate (F^) differs 



from the forward reaction rate and the net reaction rate is given by 



26 



48| 



327rG'|. cos^ ^c"^*m*/i/ 
(2^0^ 



-Rfi j deuel{J{ey - 5fi) - J{ey + 5/i)}. 



(47) 
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One important aspect of dUrca process is the opening of this channel in the forbidden 
regime Kp^ > + which was otherwise closed in field free case 48|]. The dUrca 
process also operates in the allowed domain Kp^ + Kp^ > Kp^^ in the presence of a magnetic 
field. We adopt fitting formulas for i?^ in both domains as given by Ref. 48j. 



Next we focus on the case when both protons and electrons populate zeroth Landau levels 



47|, 



48 



491]. In this case we write the matrix element as 
(27r)3 1 



W 



fi 



26 



(48) 



(49) 



In a magnetic field neutrons will be polarized because of their anomalous magnetic mo- 
ments. Hence for two different spin states of neutrons, matrix elements should be evaluated 
separately. The invariant amplitude squared is then |A^p = |A1+P + |7Vl_p, where 



ri2 



(50) 



s 

x{Wp+(A;j,)7'^(l + ^„p75)W„±(A;„)}{W;_(A;e)7M(l-7')V..(M} , (51) 



and ± signs denote the up and down spins respectively. The spinors for non-relativistic 
neutrons are given by 

where 



(52) 



(53) 



M / 

and X- = 

V V 

For non-relativistic protons in the zeroth Landau level, the spinor has the same form as 
given by Eq. fl52|) . For spin down relativistic leptons in the zeroth Landau level, the spinor 
is given by 



( \ 

1 


\ / 



(54) 
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For spin up and down neutrons, invariant amplitudes squared are 



(55) 



and 

\M-\^ = 32Gl cos^ 9,m*^m;glp{ei + pi)ie, - p,,) . (56) 

Following the same procedure as described in subsection 1111 Al and neglecting the neutrino 
momenta in momentum conserving delta functions, the final expression of forward reaction 
rate Fj is given by 

b mlm*pfii 



X 



X j deytl j dendepdeiF{en, Cp, ei)5{en - tp - ti - ty) , (57) 



where 

{\-^+\d) 5{ky),5{k^) 

Similarly we have, 

(l-^-l^).(.,),.(..) = 32G^cos^^,^' 



SG'pCos'e,{l+gnpY 1 + ^^ 1 + 



Vi 



Puz 



np 



1 + 



(58) 



(59) 



It is to be noted that z-component of neutrino momentum is smaller than its energy. We 
obtain 



F 



X 



'S{ky),S{k,) 



X J de,elJ{e,). (60) 

Now if the reverse reaction rate is FV and there is slight departure from chemical equilibrium 
then the net reaction rate is 26l |. 



e{kl - {kp^ - kpf}e-^'h 



X 



5(ky),5(k,) ' VI-' '■^^d./s(^ky),S(k,) 



j deyel{J{ey - 5fi) - J(e^ + (5/i)}. 



(61) 
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Using the following result from Ref. 26| 



de^el{J{e^ - - J(e^ + 6fi)} = — 

oU 



(62) 



we get 



G p cos Or [ 1 H 



4807r kp^kp, ^ 



-{{'^ + 9npY + alp 



X 



x{kT)^5ix. 



The zero magnetic field result is given by 



(63) 



TiiB^ = 0) 



17 
2407r 



where 



= COS^ |(1 + ^np)^ - ^) + (1 - ^"^')' - ^) - - ^i)} • 



(64) 



(65) 



IV. RESULTS AND DISCUSSION 

Nucleon-meson coupling constants of the model are obtained by reproducing the prop- 
erties of nuclear matter such as binding energy E/B = —16.3 MeV, saturation density 
uq = 0.153 fm~^, asymmetry energy coefficient aasy = 32.5 MeV and incompressibility 
K = 240 MeV and taken from Ref 52|]. The coupling strength of A hyperons with u 



mesons is determined from SU(6) symmetry of the quark model |53l. l54l. l55l|. The coupling 
strength of A hyperons to cr mesons is determined from the potential depth of A hyperons 
in normal nuclear matter 

Ua = -gaAO- + g^^AUjQ. (66) 
We take the potential depth Ua = —30 MeV as obtained from the analysis of A hypernuclei 



54 



56]. 



We adopt a profile of magnetic field given by 



B (nfe/no) = Bs + BAl-e 



(67) 



We consider different values for central field Be = 10^^ and 10^^ G whereas surface field 
strength is taken as Bg = 10^^ G in this calculation. We chose (3 = 0.01 and 7 = 3. The 



13 



magnetic field strength depends on baryon density in the above parameterization. Further 
the magnetic field at each density point is constant and uniform. The effects of anomalous 
magnetic moments of nucleons and contributions of the magnetic field to energy density 
and pressure are negligible because magnetic fields considered in this calculation are not too 
strong. 

Numbers of Landau levels populated by electrons and protons, are sensitive to the mag- 
netic field strength and baryon density. As the field strength increases, the population of 
Landau levels decreases. In a weak magnetic field, when many Landau levels are populated, 
we treat charged particles unaffected by the magnetic field. Further the effects of magnetic 
fields are most pronounced when only zeroth Landau levels are populated. Protons, elec- 
trons and muons populate zeroth Landau levels if central field strength ~ 10^^ G. Figure 
1 shows fractions of various particle species with normalised baryon density. We find large 
numbers of Landau levels of charged particles even when the magnetic field reaches its cen- 
tral value 10^^ G. Populations of charged particles are enhanced in the magnetic field due 
to Landau quantization than those of field free case (not shown in the figure). It is noted 
in Fig. 1 that the threshold density of A hyperons is shifted to I.Thq from its zero magnetic 
field value of 2.6nQ because of phase space modifications of charged particles in a magnetic 
field. 

The variation of pressure with energy density in the presence of a magnetic field with 
central field strength Be = 10^^ G (solid curve) is shown in Fig. 2. The dashed curve denotes 
the EoS without a magnetic field. The EoS in the presence of the magnetic field becomes 
stiffer when charged particles are Landau quantised. Here magnetic field contributions to 
the energy density and pressure are insignificant. 

Now we compute the relaxation time for both non-leptonic and leptonic reactions as given 



by Eq. (122|) . To calculate the matrix element we take gnp = —1.27 and gp\ = —0.72 |51| . 
and the Cabibbo angle {6c) is given by sin^c = 0.222. As we have already noted, charged 
particles populate many Landau levels in a magnetic field having central value Be = 10^^ G 
over entire density range considered in our calculation. For the non-leptonic process, when 
protons populate large number of Landau levels, we use the field free expression of F as 
given by Eq. (!36|) . For leptonic reactions we use the expression as given by Eq. (147|) when 
leptons and protons populate finite numbers of Landau levels. Chemical potentials and 
Fermi momenta of constituent particles are obtained from the EoS. The partial derivative 
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of chemical potentials with respect to baryon density can be calculated from the EoS. Using 
these inputs, we can compute relaxation times for both reactions as a function of baryon 
density at a particular temperature. Figure 3 displays relaxation time (r) of the non-leptonic 
process involving A hyperons in a magnetic field having its central value = 10^^ G and at 
different temperatures as a function of normalised baryon number density. Here r decreases 
with increasing baryon density. Further the relaxation time in a magnetic field increases 



23|. 



with decreasing temperature as was earlier noted in the field free case 

Relaxation times for dUrca reactions involving electrons and muons in a magnetic field 
with Be = 10^'' G and at different temperatures are plotted in Figs. 4 and 5 respectively. 
For leptonic processes, relaxation times are affected by the magnetic field. For the field 
free case, the dUrca process sets in at IAhq. In the magnetic field, relaxation times due to 
dUrca reactions drop sharply from large values in the forbidden domain Kp^ > Kp^, + Kp^. 
This is attributed to the behaviour of B!^^ which we discuss in details in connection with 
bulk viscosity due to dUrca processes below. The forbidden domain joins with the allowed 
domain Kp^ + Kp^ > Kp^ at a point from which relaxation times increase with baryon 
density. Like the non-leptonic case, relaxation times for dUrca processes also increase with 
decreasing temperature. 

Now we focus on the calculation of bulk viscosity due to the non-leptonic and leptonic 
processes. As soon as we know relaxation times of non-leptonic and leptonic reactions, 
we compute bulk viscosity coefficients for the respective processes from Eq. (!20|) . In this 
calculation we consider I = m = 2 r-mode and hence uo = 2/3Q. Further we take Q = 
3000s^^. In the temperature regime considered here, we have always ur << 1 for the non- 
leptonic process involving A hyperons. Therefore, we neglect that term in the denominator 
of Eq. fl20l) to calculate the hyperon bulk viscosity. The partial derivative of pressure with 
respect to neutron number density is calculated from the EoS using Eq. (123|) and the total 
derivative of neutron fraction with respect to baryon density is computed numerically from 
the EoS. As the relaxation time is a function of temperature, the bulk viscosity coefficient 
( also depends on temperature. The bulk viscosity coefficient for the non-leptonic process 
in a magnetic field with B^ = 10^^ G (dashed curve) and in the absence of a magnetic field 
(solid curve) are exhibited as a function of normalised baryon number density in Fig. 6 
at different temperatures. The non-leptonic reaction involves protons that populate many 
Landau levels in the magnetic field with B^ = 10^^ G. In this case, we adopt the field free 
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expression of the reaction rate as given by Eq. (136|) for the calculation of relaxation time 
and hyperon bulk viscosity coefficient in Eq. fl20l) . Therefore, the effects of magnetic field 
enter into hyperon bulk viscosity coefficient through the EoS which is modified by Landau 
quantization of charged particles. In Fig. 6, we find hyperon bulk viscosity in the magnetic 
field is suppressed compared with the field free case. 

We display bulk viscosity coefficient for the dUrca process in a magnetic field with Be = 
10^^ G and at a temperature T = 10^^ K as a function of normalised baryon density in Fig. 7. 
In this case electrons and protons populate many Landau levels. The dotted line represents 
the dUrca contribution in the forbidden domain Kf„ > + Kp^ ■ In this regime, reaction 
kinetics are characterised by two parameters x = ^''^ ^^'^-2%''" ^'^^ V = ^fI^^ where Np^ 
is the number of proton Landau levels. The dUrca reaction in the forbidden domain is an 
efficient process as long as a; < 10. This corresponds to baryon density < 2.3no. The large 
enhancement of bulk viscosity coefficient in this domain is attributed to the behaviour of 



p9C 



481]. It was noted i?^ = l/3atx = and it becomes very small when x > 10 [48 1. 



At a; = 0, the forbidden domain merges with the allowed domain Kp^ + Kp^ > Kp^ of the 
dUrca process. The dUrca bulk viscosity in the allowed domain is shown by the dash-dotted 
line. The result of zero field is shown by the solid line. The bulk coefficient increases with 
magnetic field in the allowed domain at higher baryon densities. 

Figure 8 and Figure 9 show bulk viscosity coefficients for dUrca processes involving elec- 
trons and muons in the presence of the magnetic field with central value = 10^^ G and at 
different temperatures as a function of normalised baryon density. In both cases contribu- 
tions to bulk viscosity coefficients due to dUrca processes come from the forbidden as well as 
allowed domains. As discussed above, the forbidden domain merges with the allowed domain 
at X = 0. For temperatures T = 10^ and 10^° K, bulk viscosity coefficients due to dUrca 
processes increase with baryon density. However the bulk viscosity for T = 10^^ K initially 
decreases and later increases with baryon density. This behaviour can be understood in the 
following way. For dUrca processes at 10^^ K, we have ujt < 1. On the other hand, we find 
ujT > 1 for dUrca processes at 10^ K and 10^° K. Consequently bulk viscosity coefficients 
have a dependence when ujt > 1 whereas it has a dependence when ujt < 1. This 
inversion of temperature dependence of dUrca bulk viscosity coefficients is not found in the 
case of hyperon bulk viscosity. 

Finally, we point out what happens in case of superstrong fields. We find that charged 
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particles populate zeroth Landau levels when Be ~ 10^^. Populations of charged particles are 
enhanced because of strong modification of their phase spaces. Further the EoS is modified 
due to magnetic field contributions to the energy density and pressure. The strong magnetic 
field enhances the hyperon bulk viscosity compared with the field free case. Similarly we note 
significant modification in bulk viscosity coefficients due to dUrca processes when leptons 
and protons are in zeroth Landau levels. However, there is no observational evidence for 
superstrong field ~ 10^^ G in neutron star's interior so far. 

V. SUMMARY 

We have investigated bulk viscosity of non-leptonic process involving A hyperons and 
dUrca processes in the presence of strong magnetic fields. In this calculation we consider 
magnetic fields with different central values = 10^^ and 10^^ G. The equation of state 
has been constructed using the relativistic field theoretical model. Many Landau levels of 
charged particles are populated for above values of central field. For a particular tempera- 
ture, the hyperon bulk viscosity coefficient is reduced compared with that of the zero field 
case. Further it is noted that the hyperon bulk viscosity decreases with increasing tem- 
perature as was earlier reported for the field free case. Bulk viscosity coefficients due to 
dUrca processes in a magnetic field have contributions from the forbidden as well as allowed 
domains. The bulk viscosity coefficients in magnetic fields having central values Be — 10^^ 
and 10^^ G are enhanced in the allowed domain at higher baryon densities than those of 
field free cases. We find an inversion of the temperature dependence of dUrca bulk viscosity 
coefficients at lO-*^^ K. We briefly discuss the effects of a superstrong magnetic held ~ lO-*^^ 
G on hyperon and dUrca bulk viscosities when zeroth Landau levels of charged particles are 
populated. However, such a superstrong magnetic field may not be a possibility in neutron 
stars. 

In this calculation, we adopt the field free hyperon bulk viscosity relation when protons 
populate large number of Landau levels. This may be an approximate treatment of the actual 
case. However the exact treatment of the effects of a magnetic field on the non-leptonic bulk 

viscosity would be worth studying when protons populate many Landau levels. Further the 
investigation of bulk viscosity in magnetic fields has important implications for the r-modes 
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in magnetars. This will be reported in a future publication. 
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Fig. 1. Fractions of different particle species in A-hyperon matter in the presence of a magnetic 
field having central value Be = 10^^ G as a function of normalised baryon density. 
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Fig. 2. Equation of state, pressure versus energy density, with a magnetic field having central 
value Be = 10^^ G (solid line) and without magnetic field (dashed curve). 
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Fig. 3. Relaxation time for the non-leptonic reaction involving A tiyperons in a magnetic field 
having central value Be = 10^^ G and at different temperatures as a function of normalised baryon 
density. 
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4. Relaxation time of dUrca reaction involving electrons in a magnetic field having central 
Be = 10^^ G and at different temperatures as a function of normalised baryon density. 
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Fig. 6. Bulk viscosity coefficient (dasfied line) for the non-leptonic processes involving A hyperons 
in a magnetic field having central value Be = 10^^ G and at different temperatures as a function 
of normalised baryon density. Field free cases are shown by solid lines. 
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Fig. 7. Bulk viscosity coefficient for tfie dUrca process involving electrons in a magnetic field 
having central value Be = 10^^ G and at a temperature 10^^ K as a function of normalised baryon 
density. The field free case is shown by the solid line. 
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Fig. 8. Bulk viscosity coefficient for tfie dUrca process involving electrons in a magnetic field 
having central value Be = 10^^ G and at different temperatures as a function of normalised baryon 
density. 
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